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Abstract 

In this paper we investigate mixing and transport in correspon- 
dence of a meandering jet. The large-scale flow field is a kinemat- 
ically assigned streamfunction. Two basic mixing mechanisms are 
considered, first separately and then combined together: determin- 
istic chaotic advection, induced by a time dependence of the flow, 
and turbulent diffusion, described by means of a stochastic model for 
particle motion. 

Rather than looking at the details of particle trajectories, fluid 
exchange is studied in terms of markovian approximations. The two- 
dimensional physical space accessible to fluid particles is subdivided 
into regions characterized by different Lagrangian behaviours. From 
the observed transitions between regions it is possible to derive a num- 
ber of relevant quantities characterizing transport and mixing in the 



studied flow regime, such as residence times, meridional mixing, cor- 
relation functions. These estimated quantities are compared with the 
corresponding ones resulting from the actual simulations. The out- 
come of the comparison suggests the possibility of describing in a 
satisfactory way at least some of the mixing properties ot the sys- 
tem through the very simplified approach of a first order markovian 
approximation, whereas other properties exhibit memory patterns of 
higher order. 

Key words: Gulf Stream, mixing, chaotic advection, turbulent diffusion, 
Markov process. 
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1 Introduction 

Western boundary current extensions typically exhibit a meandering jet-like 
now pattern. The most renowned example of this is given by the meanders of 
the Gulf Stream extension, which have been investigated in their variability 
by means of both hydrographic/currentmetric and remotely sensed data (see, 
e.g., Watts, 1983, for a survey of earlier studies; Halliwell and Mooers, f983; 
Vazquez and Watts, 1985; Cornillon et al., 1986; Tracey and Watts, 1986; 
Kontoyiannis and Watts, 1994; Lee, 1994). 

These strong currents often separate very different regions of the oceans, 
characterized by water masses which are quite different in terms of their 
physical and biogeochemical characteristics. Consequently, they are associ- 
ated with very sharp and localized property gradients; this makes the study 
of mixing processes across them particularly relevant also for interdisciplinary 
investigations. This is the case of the Gulf Stream (Bower et al., f985; Wish- 
ner and Allison, f 986; Bower and Lozier, 1994), of the Kuroshio and of the 
Brazil/Malvinas current (Backus, 1986). 

The mixing properties of passive tracers across meandering jets has been 
investigated in the recent past by a number of authors, following essentially 
two different approaches. The first one is that of dynamical models, where 
the flow is produced by integrating the equations of motion, time dependence 
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is typically produced by (barotropic or baroclinic) instability processes, and 
dissipation is present (e.g., Yang, 1996). These models account for several 
mechanisms acting in mixing in the real ocean, even if sorting out single 
processes of interest may be sometimes tricky. 

A second and simpler approach, the one followed in this paper, is that of 
kinematic models (Bower, 1991; Samelson, 1992; Dutkiewicz et al., 1993 - 
hereafter respectively referred to as B91, S92, DG093 - Duan and Wiggins, 
1996; for slightly different kinds of flows see also Lacorata et al., 1995). 
In such models the large-scale velocity field is represented by an assigned 
flow whose spatial and temporal characteristics mimic those observed in the 
ocean. However, the flow field may not be dynamically consistent in the sense 
of being a solution of the equations of motion, or of conserving, e.g., potential 
vorticity. Despite their somehow artificial character, these simplified models 
enable to focus on very basic mixing mechanisms. 

The paper B91 represents a first attempt at understanding particle ex- 
change in a two-dimensional meandering jet steadily propagating eastward. 

The large-scale flow proposed in B91 has been utilized as a background 
field in further works where mixing is separately enhanced by two different 
transport mechanisms. S92 considers a modification of the B91 flow field 
where fluid exchange is induced by chaotic advection generated by a flow 
time dependence. The basic flow is made time-dependent in three different 
fashions: the superposition of a time-dependent meridional velocity; that of a 
propagating plane wave; a time oscillation of the meander amplitude, which 
is the case we further investigate in this paper. 

The Melnikov method (see Lichtenberg and Lieberman 1992, LL92 here- 
after) is used in S92 to explore the chaotic behaviour around the separatrices 
of the original B91 flow when a time dependence is added. One of the results 
of this investigation is that while mixing occurs between adjacent regions, 
over a broad range of the meander oscillation frequencies, it does not easily 
take place across the jet, i.e. from recirculations south of the jet to recircu- 
lations north of it. This is inherently due to the oscillation pattern of the 
large-scale velocity, and we will discuss this in detail further in this paper. 

Particle exchange in the same B91 flow is achieved by DG093 by super- 
imposing to the original time-independent basic flow a stochastic term which 
describes mesoscale turbulent diffusion in the upper ocean. The focus of that 
paper is on the exchange among recirculations and jet core and viceversa, and 
on the homogenization processes in the recirculation. The numerical experi- 
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ments presented in DG093 are carried out for quite short integration times, 
which do not allow for exploring the mixing across the jet. 

Since in the real ocean the two above mixing mechanisms, i.e. chaotic 
advection and turbulent diffusion, are simultaneously present, in this paper 
we investigate how particle exchange varies through the progression from pe- 
riodic to stochastic disturbances, revisiting and putting together the mixing 
processes studied by S92 and DG093. 

This is done by looking at particle statistics obtained by numerical com- 
putation of the trajectories of a large number of particles (or equivalently, 
since our system is ergodic, following one particle for a very long time) in 
three different flows: one equivalent to S92, in which mixing is induced by 
chaotic advection; one equivalent to DG093, where it is due to turbulent 
diffusion, and a combination of them. 

Dispersion processes in a flow field can be quantitatively characterized, 
in the Lagrangian description, in terms of different quantities, such as, e.g., 
the Lyapunov exponent A (Benettin et al. 1980) and the diffusion coefficients 
Dij (LL92). 

However the above indicators, even if mathematically well defined, can 
be rather unrelevant for many purposes. The Lyapunov exponent is the in- 
verse of a characteristic time t^, related to the exponential growth of the 
distance between two trajectories initially very close; however, other charac- 
teristic time scales may appear and result as relevant in the description of a 
system, such as those involved in the correlation functions and in the mixing 
phenomena. It is worth stressing that there is not a clear relationship, if any, 
among these times and t^. 

Also the use of the diffusion coefficients can have severe limitations; some- 
times the are not able to take into account the basic mechanisms of the 
spreading and mixing (Artale et al., 1997). Our western boundary current 
extension system has essentially a periodic structure in the zonal direction 
It is thus possible to define and (numerically) compute the diffusion coeffi- 
cients. They are related to the asymptotic behaviour, i.e. long times and 
large spatial scales, of a cloud of test particles. On the other hand, if one is 
interested in the meridional mixing, which takes place over finite time scales, 
the diffusion coefficients may not be very useful. In such situations it is then 
worthwhile to look for alternative methods of describing mixing processes, 
as was done by Artale et al. (1997) looking at dispersion in closed basins or 
by Buffoni et al. (1997), who employ exit times for the characterization of 
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transport in basins with complicated geometry. 

Our investigation is carried out with a non conventional approach, in a 
geophysical contest, as we try to analyze the system from the standpoint of 
the approximation in terms of markovian processes (Cecconi and Vulpiani 
1995, Kluiving et al 1992, Nicolis et al. 1997, Fraedrich and Muller 1983 and 
Fraedrich 1988). 

We start from the consideration that the flow field we want to character- 
ize in terms of fluid transport can be subdivided into regions corresponding 
to different lagrangian behaviours: ballistic fly in the meandering jet core, 
trapping inside recirculations, retrograde motion in the far field. As an ob- 
vious consequence, we introduce a partition of the two-dimensional physical 
space accessible to fluid particles and divide it into these disjoint regions se- 
lected in a natural way by the dynamics. At this point, one can stydy the 
transition of fluid particles between different regions as a discrete stochastic 
process generated by the dynamics itself. 

In this paper we study the statistical properties of this stochastic process 
and we compare it with an approximation in terms of Markov chains. For 
some fluid exchange properties — such as, e.g., the probability distribution 
of the particle exit times from the jet or from the neigbouring recirculation 
regions — the effects of the two different mixing mechanisms and the results 
of the markovian approximation are very similar. Other properties, such as 
the meridional mixing across the jet, do not show such an obvious possibility 
to be described in terms of markovian simple processes. However, for those 
properties the markovian description is seen to be relatively more accurate 
in the case when chaotic advection and turbulent diffusion are simultaneosly 
present. 

The comparison between the results of the numerical simulations and 
those computed in the markovian approximation allows for a deeper under- 
standing of the transport and mixing mechanisms. 

In section 2 we introduce the kinematic model for the flow field corre- 
spondent to the Gulf Stream flow and both models for chaotic advection and 
turbulent diffusion. Section 3 is devoted to the description of the markovian 
approximation. In Section 4 we discuss the numerical results and the com- 
parison of the true dynamics with the markovian approximation. Section 5 
contains some discussion and conclusions. The Appendix summarizes some 
basic properties of Markov chains. 
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2 The flow field 



The large-scale flow in its basic form, representing the velocity field in cor- 
respondence of a meandering jet, is the same introduced in B91 and further 
discussed in S92; in a fixed reference frame the streamfunction is given by 
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where ipQ represents half of the total transport, A, k and c x the amplitude, 
wavenumber and phase speed of the pattern. A change of coordinates into 
a reference frame moving eastward with a velocity coinciding with the phase 
speed c x , and a successive nondimensionalization, yield a streamfunction as 
follows: 
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The evolution of the tracer particles is given by: 



dx 
dt 



dy 



dy 
dt 



dx 



(3) 



The natural distance unit for our system is given by the jetwidth A, set 
to 40 Km (B91, S92). The basic flow configuration is very similar to case (b) 
of B91: B was chosen as 1.2; c as 0.12. The only sensible difference is the 
value assigned to L, i.e. the meander wavelength, which was set as 7.5, as 
will be discussed below. 

In fig. 1 we show the stationary velocity vector field in the moving frame: 
the field is evidently divided into three very different flow regions (see also 
B91, S92, DG093): the central, eastward moving jet stream, recirculation 
regions north and south of it, and a far field; the far field, given our choice of 
parameters, appears to be moving westward at a phase speed of — c x = —0.12. 
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This intrinsic self subdivision of the flow field will result crucial for building 
a partition of the possible states available to our test particles, which will be 
investigated in markovian terms. 

Chaotic advection may be induced in a two-dimensional flow field by 
introducing a time dependence (see e.g. Crisanti et al., 1991). This is simply 
achieved by adding to the basic steady flow some tipically small perturbation 
which varies in time. Among three basic mechanisms discussed by S92, we 
chose here a time-dependent oscillation of the meander amplitude: 

B(t) = B + e • cos{ut + 9) (4) 

In (|j) we set B = 1.2, e = 0.3, uj = 0.4 and 9 = tt/2. These choices, as 
well as that for L, are motivated mainly by the results of the observations by 
Kontoyiannis and Watts (1994) and of the numerical simulations by Dimas 
and Triantafylou (1995). Namely, the most unstable waves produced in the 
latter work compare very well with the observations of the former, which show 
wavelengths of 260 Km, periods of ~ 8 days, e-folding space and time scales 
of 250 Km and 3 days respectively. In our case, since the downstream speed 
was set to lm/s, our e-folding time scale would correspond, in dimensional 
units, to approximately 3.5 days. The flow field resulting from the time 
dependent version of (0) is shown in fig. 2 for three different subsequent 
time snapshots t = 0. (T/2) (fig. 2a), t = T/4 (fig. 2b) and t = 3T/4 
(fig. 2c). Our system shows two different separatrices with a spatial periodic 
structure (see fig. 2), North and South of the jet. At small e one has chaotic 
motion around them but without meridional mixing. In order to have a 
"large scale chaos" , i.e. the possibility that a test particle passes from North 
to South (and viceversa), crossing the jet, one needs the overlapping of the 
resonances (Chirikov 1979) e > e c . In our case e > e c and e c depends on u 
(in fig. 3 we show e c vs uo for our system). 

The lagrangian motion of a test particle is formally a Hamiltonian system 
whose Hamiltonian is the stream function 0. If = <p (x,y) + 5<p(x,y,t), 
where S(p(x, y, t) = 0(e) is a periodic function of t, there exists a well known 
technique, due to Poincare and Melnikov, which allows to prove whether the 
motion is chaotic (LL92). Basically if the steady part <f)o(x,y) of the stream 
functions admits homoclinic (or heteroclinic) orbits, i.e. separatrices, then 
the motion is usually chaotic in a small region around the separatrices for 
small values of e. 
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In S92 the Melnikov integral has been computed explicitly for the 0o 
and 5(fi that we have used in this paper, and thus proved the existence of 
lagrangian chaos. However, even if the Melnikov method can determine if 
the lagrangian motion is chaotic, it is not powerful enough for the study of 
other interesting properties which will be the focus of section 4.1. 

Alternatively (or jointly), mixing in the flow field (|2|) can be created by 
adding a turbulent diffusion term. This was done utilizing a stochastic model 
for particle motion belonging to the category of the so-called "random flight" 
models (e.g. Thomson, 1987), which can be seen as simple examples of a more 
general class of stochastic models which can be nonlinear and have arbitrary 
dimensions, described by the generalized Langevin equations (Risken, 1989; 
for a review, see Pope, 1994): 

dsi = hi(s,t)dt + g i j(s,t)d[i j [i = l,...,N] (5) 

where s = (si, . . . , sn) are N stochastic variables which, in our context, are 
the turbulent velocity fluctuations, /i, is a random process with independent 
increments, and hi and g^j are continuous functions. A general, remarkable 
characteristic of these models is their markovian nature, which obviously has 
a particular interest for this investigation. The theoretical motivation for the 
choice of markovian models to describe mesoscale ocean turbulence has been 
thoroughly discussed in Zambianchi and Griffa (1994a), Griffa (1996) and La- 
corata et al. (1996); it is worth adding that this particular model has proved 
to accurately represent upper ocean turbulence in regions characterized by 
homogeneity and stationarity (see Zambianchi and Griffa, 1994b; Griffa et 
al., 1995; Bauer et al., 1997; but also the results of numerical simulations by 
Verron and Nguyen, 1989; Yeung and Pope, 1989; Davis, 1991), and is easily 
extended to more complex situations (van Dop et al., 1985; Thomson, 1986). 

In our simulations, a turbulent velocity <5u^(x, t) is added to the large- 
scale velocity field u^(x, t) resulting from the streamfunction (0). The 
resulting equation for the particle trajectory is: 

§ = u(x, t ) (6) 

where u(x, t) is given by 

u(x, t) = u (M) (x, t) + 5u (T) (x, t). (7) 
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Our model assumes 5vS T ) (x, t) as a Gaussian process with zero mean and 
correlation: 

< £«f ) (x,t)5ufV,0 >= 2a% j 5(x - ^e^'^. (8) 

With this choice, 5u^(x, t) is a linear in time markovian process: the tur- 
bulent field is entirely described in terms of two parameters: the variance of 
the small scale velocity a 2 and the e-folding time scale of the velocity auto- 
correlation function, i.e. its typical correlation time scale r, which represents 
the time step of the markovian process. The interdependence among smaller 
and larger time scales of the lagrangian motion will be investigated in the 
following chapters. 



3 The Markovian approximation 

The idea to use stochastic processes to study (and describe) chaotic behaviour 
is rather old (Chirikov 1979, Benettin 1984). One of the most relevant and 
successful approach is the symbolic dynamic, which allows to give a detailed 
description of the statistical properties of a chaotic system in terms of a 
suitable discrete stochastic process (Beck and Schogl 1993). 
Given a discrete dynamical system: 

x t = S'x , (9) 

one can introduce a partition A dividing the phase space in Ai, A 2 , • ■ ■ , A^ 
disjoint sets (with Ai fl A,- = if i ^ j). From each initial condition one has 
a trajectory: 

x ,xi,---,x n , ••• (10) 

The point x G A io will select the integer i , the next one Xi G A ix the integer 
ii and so on. Therefore for any initial condition x we have a certain symbol 
sequence: 

x -»• (i ,ii, •••»»!»,•••)■ (11) 
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Now the study of the coarse grained properties of the chaotic trajecto- 
ries is reduced to the statistical features of the discrete stochastic process 
(i , i%, ■ • • , i n , • • •)■ A useful and important characterization of the properties 
of symbolic sequences is the Kolmogorov- Sinai entropy, defined by: 

h KS = \iva{H n+1 -H n ), (12) 

n— »oo 

with 

H n = sup[- £ P(Cn) In P(C n )], (13) 
Ml c n 

and 

C n = (io,h, ■ ■ ■ ,i n -i) , (14) 

where P(C n ) is the probability of the sequence C n and {.4} is the set of all 
possible partitions. 

We note that, from a theoretical point of view, the sup in (|T3| ) hides a very 
subtle point: sometimes there exist a particular partition, called generating 
partition, for which one has automatically the sup. A partition is generating 
if the infinite symbol sequence iq, i±, • • • , i n , • • • uniquely determines the initial 
value xo. Unfortunately it is not trivial at all to know if a system possesses a 
generating partition. Moreover in practical applications, even if the system 
admits a generating partition, it is extremely hard to find it (see Beck and 
Schogl 1993 for more detail). 

A part the difficulties for the choice of a suitable partition, the stochas- 
tic process given by the symbol dynamics with a given partition can have 
rather nontrivial features. Of course the optimal case is when the symbolic 
stochastic process is a Markov chain, i.e. the probability to be in the cell A{ 
at time t depend only on the cell at time t — 1. In this case is possible to 
derive all the statistical properties (e.g. entropy and correlation functions) 
from the transition matrix Wy whose elements are the probabilities to have 
the system in the cell Aj at time t if at time t — 1 the system is in the cell 
Aj. See the Appendix for a summary of the properties of Markov chains. 

Usually a /c-order markovian process (i.e. one in which the probability 
to be in the cell Aj at time t depends only on the preceding k steps t — 
1, t — 2, ■ • • , t — k), with large k, is necessary to mimic, with a good accuracy, 
a chaotic system. In particular by means of the quantities defined in eqs. 
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( |T2| , P~3D and of consideration of information theory (see Khinchin 1957) we 
are able to estimate the order of the Markov process necessary to reproduce 
the statistics of the process (i , i±, • • • , i n , ■ • ■) generated by the dynamics. In 
fact it is possible to demonstrate (Khinchin 1957) that, defining 

h n = H n -H n _ 1 (15) 

with H n _i given by flT3|), if the process is a Markov process of order k then 
h n = hxs for each n > k + 1. In the next section we apply this method, in 
our case, to give an estimate of the order k. 

However we observe that usually it is relatively easy, using a markovian 
approximation of order k < 4 4-5, to find a reasonable agreement for the 
K-S entropy ([12]) (or for the Lyapunov exponent). On the other hand for 



different properties, such as the correlation functions, it is necessary to use 
Markov processes of rather large order (k ~ 8 4- 10) just for a fair agreement. 

The issue to mimic a low dimensional dynamical system in terms of a 
Markov process of a certain order is surely an interesting aspect in the field 
of chaotic dynamics but, in our opinion, with a weak relevance for many 
practical purposes in geophysics since one needs a very large statistics for 
the computation of the transition probabilities. Therefore we shall restrict 
our analysis to the simplest case of the approximation in terms of a Markov 
chain, i.e. a first order process. This practical approach has been successfully 
used in the study of certain problems related to the dynamical properties of 
small astronomical bodies such as comets (see Rickmann and Froeshle, 1979 
and Levinson 1991) and for the interpretation of geophysical phenomena 
(Fraedrich, 1988, Fraedrich and Miiller, 1983). 

We shall now explain how we proceeded to express the behaviour of our 
Gulf Stream-like system in terms of symbolic cynamics. First we reduced 
the ordinary differential equation (|^) obtained by the stream function (Q) 
to a discrete in time dynamical system, which was accomplished building 
the Poincare map associated with (||]). This represents an assessment of the 
possibility to write x n+1 in terms of x n - defined as x n = x(t =nT) -: 



F[x n ] . (16) 



Writing down an explicit expression for F[x n ] may be in general non trivial; 
however, it is worth noticing the importance the existence of such a relation- 
ship. 
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We have now to decide a suitable choice of the phase space partition. Dif- 
ferently from the case of the orbital evolution of comets, where the cells are 
rectangular regions of the phase space, for our system we shall adopt a parti- 
tion whose cells have curved frontiers. Considering the streamline pattern of 
our flow field (fig. 2), the structure of the physical space accessible to fluid 
particles suggests an obvious, natural choice for the partition: a particle will 
find itself in state 1 when it is inside the jet core (open trajectories); states 
2 and 3 correspond to trapping in the Northern and Southern recirculations 
respectively (closed trajectories) and states 4 and 5 to the far field, i.e. far 
from the jet (backward open trajectories). 

This partition turns out to be particularly appropriate to describe some 
important mixing properties in our system, such as: 

• (a) The residence times of particles in the trapping recirculations or 
inside the jet, which in the language of Markov chain (see Appendix) 
correspond to the first exit times from state i {i = 1, • • • , 5); 

• (b) The meridional mixing times (MMT), i.e. the time it takes to a 
particle to enter the Northern (Southern) recirculation starting from 
the Southern (Northern) one, i.e. the time of first passage from cell 2 
(3) to cell 3 (2); 

• (c) The correlation function for a variable Xi( n ) which indicates if a 



determined state i is visited at time n (see below, eqs ^3] and ^4], and 
Appendix). 

Because of the system symmetries, states 2 and 3 possess the same statistical 
properties and so do states 4 and 5; in particular, the following equalities 
hold: 

W12 = Wi 3 , W 23 = W 32 , W 21 = W 31 , W 22 = W 33 , and so on (17) 

Let us now describe how to compute statistics for the quantities (a,b,c). 
First of all, we can compute from a long trajectory xo,xi, • • • ,x n (n ^> 1) 
the transition probabilities: 

W„ = lim %M (18 ) 

n-^oo N n (l) 
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where N n (i) is the number of times that, along the trajectory, Xf (t < n) visits 
the cell Ai and N n (i,j) is the number of times that x 4 e Aj and x t+ i e A;. 
In table 1 we report the elements of the matrix. 
Notice that, for each % : 

E^ = 1 (19) 

3 

We can express the probabilities {Pi} in terms of the matrix {Wjj}: 

Pi = J2 P 3 W ji (2°) 

3 

Let us stress that the eqs. (|I5| , p0| ) hold even if the process is not a Markov 
chain. Under the assumption (approximation) that the symbolic stochastic 
process generated by our deterministic chaotic model is a Markov chain one 
can derive (see Appendix) the probability of the first exit times from state % 
in n steps: 

" i • - uv: 



Pn) 



W, 



mr, (2i) 



which is the statistics of residence times in state i. A slightly more compli- 
cated computation gives the probabilities fij{n) of first passage from state i 
to state j in n steps: 



n-l 



f l3 {n) = {W n ) t3 - J2 fij(n - k)(W% , (22) 

k=l 

i.e. the statistics of the MMT. For the normalized correlation function Cj(ra) 
of the variable Xi{ n ) defined as: 

X ^ = { 'otherwise (23) 

we have: 

C, l (n) = ((W n ) ll -P l )/(l-P l ). (24) 

The Kolmogorov- Sinai entropy for the Markov chain is nothing but the 
Shannon entropy for a Markov chain: 

h K s = h s = -J2 P i W H ln W ir (25) 
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Note that for a Markov process h n = H n — H n -i = hxs for n = 2 (see above). 
Since the discrete time system is obtained looking it at times 0, T, 2T, ■ ■ ■ the 
Lyapunov exponent A of the original system has to be compared with: 

Am = J,' ( 26 ) 

This last equation is easily understood noticing that hxs gives the degree 
of information per step produced by the process, that for a chaotic system 
in two dimensions corresponds to the Lyapunov exponent apart from a time 
rescaling. 



4 Numerical results 

We now discuss the numerical results for the models introduced in section 2 
and their comparison with the markovian approximation illustrated in section 
3. 



4.1 Mixing induced by chaotic advection 

We first consider the deterministic model with the parameter B of the stream 
function eq. @ varying periodically in time according to eq. @ with the 
parameters B = 1.2, e = 0.3, u = 0.4, <p = tt/2 and c = 0.12. With this 
choice the system is chaotic and exhibits mixing at large scale, i.e. North- 
South mixing occurs. 

We show in fig. 4 the spreading at different times of a cloud of test 
particles. The domain is naturally defined from the basic cell that repeats 
itself creating a zonal periodic structure of wavelength L; x thus varies in 
[0, L\ , while y in [—4, 4] . We fixed a posteriori these bounds for y since for 
our choice of parameters no particles reach the far field, and no trajectories 
attains values in \y\ larger than |4| (even though, in general, we expect low 
but non zero frequencies for these states, see also S92). 

In general, whether a North-South mixing happens or not depends sensi- 
bly on the values of e and cu. Typically the system reveals a strong preference 
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to have a long residence times in the Northern or the Southern half of the 
domain with respect to the jet core. This peculiar feature will play an im- 
portant role in the comparison with the markovian approximation. 

The transition matrix elements Wij and the visit probabilities {Pi} are 
computed by means of eq. fll8|) , looking at x every period i.e. for t = 
T, 2T, • ■ ■ where T = 2tt/uj, see table 1 and 2. At a first glance, we can see 
that the requested symmetry properties are respected (eq. |17|) . 

Now in order to test whether the system is well approximated by a first 
order Markov process we compute exit times, correlation functions, merid- 
ional mixing times and Lyapunov exponent from the actual dynamics and 
compare them with the markovian predictions. 

In fig. 5 we show the first exit time probability distributions for the states 
1, 2 (3) and the corresponding markovian predictions. After stressing that the 
straight lines of fig. 5 are not to be confused with best-fit curves, we see that 
the agreement is good over a certain range both for state 1 and 2 (3). The 
agreement between the markovian predictions and the actual results is rather 
poor for small and very large exit times. The above behaviour shows that 
the markovian approximation cannot hold at small times since the details 
of the dynamics are strongly relevant. In the same way non trivial long 
time correlations cannot be accounted for in terms of a first order markovian 
process. 

In fig. 6 we can see how the correlation functions of Xi ( see ec L s - p3[-p4|) 
for the state 1 and 2 (3) are just in vague agreement with the corresponding 
correlations obtained from the markovian process. The trajectories in the 
recirculations (i.e. states 2 and 3) appear to be much more auto-correlated 
than those in the jet. Therefore we deduce that the system, although chaotic, 
has a strong memory as to which half (North or South) of the spatial domain 
it is visiting. So the typical evolution is a "rebound game" between state 2 
and the South half of the jet during a certain time interval; then it crosses 
the jet core and performs again the same pattern between state 3 and the 
North half of the jet, until it jumps back; and so on. 

This is strikingly evident looking at the distribution probabilities of the 
meridional mixing times, in fact in this case the markovian approximation 
completely fails (fig. 7). 

The markovian approximation leads to a clear disagreement in this case 
because it actually happens that when a tracer leaves a recirculation region, 
say in the South half, and goes into the jet, most of the time it returns back to 
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some other Southern close orbit rather than passing through the jet barrier, 
while in the first order Markov approximation it is not possible to explain 
such a strong memory effect. This feature is a clear indication that higher 
order Markov processes are necessary to describe the statistics produced by 
the dynamics. This is shown in fig. 8, where transitions between states 1 
and 2, and 1 and 3, are compared: whereas in the first order markovian 
approximation (solid lines) a test particle jumps very often from North to 
South, the results of our numerical experiments show a stronger tendency for 
particles to keep being confined either between states 1 and 2 or between 1 
and 3. 

In order to quantify the relevance of the memory effects we computed the 
block entropies h n , defined in section 3 (eq. |T5|) at varying n. In fig. 9 we 
can see that to obtain the convergence of the entropies we need at least of a 
Markov approximation of order 6 -v- 7. 

The Lyapunov exponent computed with a standard algorithm (see Benet- 
tin et al. 1980) is A = 0.05, the first order Markov approximation gives 
Aju = 0.03, while the extrapolation with the asymptotic value h = lim^oo h n 
gives Am = 0.03. The fact that Am < A is probably due to the fact that the 
partition here used is not a generating one (see section 3) however there exist 
a fair agreement between Am and A. It is worth noticing that the above fea- 
tures are fairly robust, and do not vary in a relevant way after weak changes 
of parameters. 



4.2 Mixing induced by turbulent diffusion 

Now we perform the same investigation discussed in the previous subsection 
setting e = and turning on turbulent diffusion which is described in terms 
of a stochastic model for particle motion: 

dx dy , s 

H = u + Vl > Tt =v + T]2 > (27) 

where u, v are given by the stream function (|^) and rji are zero mean gaussian 
stochastic processes with < r]i{t)r]j{t') >= a 2 5ij exp[— \t— t'\/r]. The gaussian 
variable rji are generated by a Langevin equation (Chandrasekhar, 1943): 



drjj 
dt 




1,2 (28) 
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where the variables Q are zero mean Gaussian noises with < d{t)(j (£') >= 
5ij5(t — t'). The numerical integration of the equations has been performed 
by means of a stochastic 4-th order Runge-Kutta algorithm (Mannella et al., 
1989). 

Now the motion is unbounded for the presence of the isotropic diffusive 
terms, so we need a "trick" to limit the dispersion along y inside a domain 
as similar as possible to the previous one. In this way we can hope to look 
for a transition matrix comparable with that one obtained in the chaotic 
deterministic case and compare the two models. 

Since the chaotic model does not fill the states 4 and 5, we impose that 
if a tracer enters a backward motion region it is reflected back by changing 
the sign of the meridional turbulent velocity, being the stream lines in these 
regions almost parallel to the x direction. Let us stress that this boundary 
condition practically does not affect the mixing processes between the jet 
and the recirculation regions, therefore we simply need to follow only those 
branches of trajectory which fall within the domain of study and not the 
large scale diffusion motion far from the stream frame. 

In the diffusive case we fix the values a = 0.05 and r = T/4 as rep- 
resentative of an observable situation (see DG093) and compute again the 
transition matrix and stationary frequencies of the 5 states (see tables 1 and 
2). 

In fig. 10 we show the spreading of a cloud of test particles. Notice that 
the artifact introduced to bound the dynamics in the y direction does not 
produce any artificial accumulation of particles near the boundary. 

At first we notice that the elements of the transition matrix are close to 
those of the chaotic case. We have computed the transition probabilities over 
a time period T as for the chaotic case so that we can reasonably compare 
the two cases. 

Fig. 11 shows the probabilities of the first exit times of the states 1, 2 
and 3 and the relative markovian predictions. These distributions are very 
well approximated by the first order Markov process. 

The correlation functions are shown in Fig. 12. The difference between 
the actual and the markovian curves is now smaller than in the chaotic case 
because of the presence of diffusion that decreases sensibly the degree of 
memory. 

The most relevant difference with respect to the chaotic model is a clear 
improvement in the agreement with the markovian approximation for the 
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distributions of the meridional mixing times, see fig. 13. The memory effects 
are now smoothened by diffusion, and the transition rates are much more 
representative of the actual dynamics. Thus we can conclude that in the 
diffusive model the markovian approximation works much better than in the 
chaotic one. Looking at the block entropies h n (|l^) we have a clear indication 
that the process is of a lower order respect the chaotic case (compare fig. 14 
with fig. 9) 

Just like in the previous case we have investigated the behaviour of the 
system varying the parameters a and r. We have observed that if we keep 
the quantity a 2 r constant the system displays similar behaviours, even if 
the extent of the agreement between simulation results and markovian ap- 
proximation slightly differs for different values of the turbulence parameters; 
this can be understood if we recognize that a 2 r corresponds to the diffusion 
coefficient (see, e.g., Zambianchi and Griffa, 1994a). It has been shown that 
varying a and r even if keeping the diffusion coefficient constant can in- 
deed affect the quantitative estimates of dispersion in cases characterized by 
inhomogeneity and/or nonstationariety (see, again, Zambianchi and Griffa, 
1994a). On the other hand, the qualitative functional behaviour of the dis- 
persion processes has been seen to be influenced very little by such changes 
in the parameters of turbulence (Lacorata et al., 1996). 



4.3 Mixing jointly induced by chaotic advection and 

turbulent diffusion 

A detailed analysis of Lagrangian data from the ocean aimed at determin- 
ing contemporary presence and relative importance of chaotic and turbulent 
mixing is at present still lacking, as it would imply the evaluation of both 
one- and two-particle statistics parameters, which has been done only in the 
context of purely diffusive particle exchange investigations (see, e.g., Poulain 
and Niiler, 1989). However, in the real ocean we expect both the above 
mixing mechanisms, discussed in sect. 2, to be present at the same time. 

One of the interesting results of the previous sections is not only the 
fact that one can look at particle exchange in terms of markovian processes, 
but also that the sampling time suitable for the description of chaotic and 
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turbulent exchange are of the same order of magnitude for fairly realistic 
simulations. This suggests the feasibility of a numerical experiment in which 
a stochastic term is added to a time dependent large-scale velocity field. 

In addition, in the Introduction we mentioned the issue of a possible in- 
consistency of kinematic models as to the lack of lagrangian conservation of 
quantities such as potential vorticity. This difficulty, which has been dis- 
cussed at length recently for two-dimensional chaotic flows (see, above all, 
Brown and Samelson, 1994; Balasuriya and Jones, 1996), is overcome in the 
combination of the two mixing processes, as turbulent diffusion can be seen as 
a sort of dissipation, which therefore acts so as to "smear" potential vorticity 
gradients. 

In this numerical experiment we use the model equations (|5|, |7]) with 
U ( M ) given by the stream function (0) with the time- dependent pertur- 
bation (||) and for the turbulent velocity (W T ) (x, t) we use the stochastic 
process defined in eqs. ( |2T| - [28D ; the parameters are: B = 1.2, e = 0.3, uj = 
0.4, 9 = 7r/2 and a = 0.05, r = T/4 where T = 2ti/lj. As one can see this 
choice for the parameters is simply a superposition of the two previous cases. 

Also in this case the matrix elements (i.e. the transition probabilities) are 
comparable we the other ones (see Table 1). As we can see from figs. 15 and 
16 the distributions of the residence times and of the meridional mixing times 
display the same qualitative behaviour as the pure diffusive case (compare 
with figs. 11 and 13), from which we can deduce that for these features the 
most relevant effect is due to the diffusive term. 

As to the correlation function (fig. 17) we note that there is a sensible 
improvement for the markovian approximation with respect to either the 
purely chaotic or the purely diffusive case, since the superposition of the 
two independent perturbations to the large scale flow strongly decreases the 
memory effects. 



5 Discussion and conclusions 

In this paper particle exchange in a meandering jet has been investigated 
by means of a kinematic model in which mixing is obtained by two differ- 
ent mechanisms: chaotic advection and turbulent diffusion. The large-scale 
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structure of the jet-like flow is assigned in terms of a stationary streamfunc- 
tion. This has been modified in two ways, in order to provide the requested 
fluid exchange: chaotic advection is induced by adding a time-dependent, 
relatively small perturbation to the steady portion of the streamfunction. 
Alternatively, turbulent diffusion has been introduced by superimposing to 
this latter a stochastic field. The turbulent field has been selected so as to 
resemble as closely as possible the typical effect of upper ocean turbulence in 
the absence of coherent structures. Numerical simulations have been carried 
out for a case in which the two above effects have been jointly present, trying 
to take into account the richness and complexity of situations observed in the 
ocean, where the two different mixing mechanisms are thought to be present 
simultaneously, even if possibly acting at different time and space scales. 

The instrinsically different nature of the two investigated mixing mech- 
anisms has resulted in the past in disjointed descriptions of their respective 
effects: chaotic advection in correspondence in meandering jets has been 
studied, e.g., by means of methods derived from the dynamical systems the- 
ory (Pierrehumbert, 1991; S92; Wiggins, 1992; Duan and Wiggins, 1996), 
whereas the action of turbulent diffusion was clarified by phenomenological 
Lagrangian motion analysis (B91, DG093). 

In this paper mixing is studied in terms of particle transitions among areas 
of the physical two-dimensional space characterized by qualitatively differ- 
ent flow regimes, observed as realizations of a markovian process. Given the 
structure of the velocity field, the partition of the space accessible to parti- 
cles is self-evident and physically consistent. A delicate point is obviously the 
choice of the appropriate scale for the sampling time of the process. However, 
in our cases an inherent time scale is present in the velocity field, and this is 
set by, in turn, either the space-time structure of the deterministic portion of 
the flow (chaotic case) or by the intrinsic memory time scale present in the 
stochastic velocity. The markovian approach results, in this sense, a quite 
natural one to undertake looking at the overall mixing from a unified per- 
spective, embedding elements of dynamical systems and of stochastic process 
theory. Also, it is an alternative way to look at diffusion avoiding the usual 
diffusion coefficients, whose general relevance in geophysics has been recently 
subject to debate (see Artale et al.1997 and Buffoni et al. 1997). 

For some fluid exchange properties, the effects of the two above mixing 
mechanisms are comparable with the results of the markovian approxima- 
tion: this is the case, for instance, of the exit times of particles from the 
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jet and the recirculating regions north and south of it. On the other hand, 
chaotic advection and turbulent diffusion act quite differently, under that 
perspective, when it comes to meridional mixing and correlation functions. 
The failure of the markovian approximation for the characterization of the 
meridional particle exchange in the chaotic is due to non-trivial long term 
memory effects. Since turbulent diffusion is modelled by a non-white noise 
process in the stochastic velocity field, we would expect for the turbulent 
case a closer behaviour to that predicted by the markovian approximation. 
For the same reason, given the results for the purely chaotic simulations, the 
combined effect of chaos and diffusion has been expected to be well described 
in markovian terms. This is indeed the case, and the results for the joint fluid 
exchange situations agree quite closely with the markovian predictions. 

It is worth underlining that the possibility to look in terms of a markovian 
approximation at mixing in regions characterized by a quite complex flow 
structure, even in the presence of different transport mechanisms, can have 
quite interesting applicative consequences. When the small-scale details of 
mixing are beyond our interest, and if and when our flow system shows fairly 
well defined time scales, it is apparently possible to look at particle exchange 
in a relatively simple manner, over time scales which allow for a sensible 
reduction of the sampling rate. This aspect often turns out to be a critical 
constraint for the undertaking, e.g., of Lagrangian investigations of the real 
ocean, where reducing the required sampling rate can result in reducing the 
amount of data to be collected and transmitted. 
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Appendix 



A Some properties of Markov chains 

A Markov chain is a stochastic process in which the random variable de- 
scribing the state of the system (in our case the cell occupied by the tracer) 
and the time are discrete and the probability to be in a given state at time 
n depends only on the state at time n — 1. All the properties of a Markov 
chain can be derived from the transition matrix, {W 7 ^}, whose elements are 
the probabilities to be in state j at some time n being at time n — 1 in state 
i. For example the probability to go to state j starting from i in n steps is 
simply: 

Prob(i -> j;n) = {W n )ij . (29) 

An excellent introduction to Markov chain can be found in Feller (1968). 
In this appendix we only resume some formulas that are useful in describing 
some relevant properties of our system. 

First of all we have: 

E w« = 1 • (3°) 

3 

In addition to matrix {Wy} one can compute the stationary probabilities Pj 
to visit the cells Ai as elements of the (left) eigenvector corresponding to the 
eigenvalue 1 : 

P^J^PiWij, (31) 

i 

Let us note that eqs. (^,^T|) are rather general results that hold for a generic 
discrete stochastic process. Eq. (|30| ) describes the conservation of probability 



and eq. (plf ) is nothing but Bayes' Theorem. In order to have ergodicity and 
mixing properties the Markov chain must have a non zero probability to pass 
through any state in a finite number of steps, i.e. there exist a n such that 
(W n )ij > (Feller 1968). Defining pi(t) the probabilities to visit the state i 
at time t, for a Markov chain we have: 

Pi(t+l) = X>(*)W« (32) 
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in this way we have the evolution of the probability vector pi (see Rickman 
and Froeschle, 1979), eq. flST]) corresponds to t — > oo of eq. (|32D, i.e. the 
equilibrium distribution pj(oo) = P { . 

The probability of the first exit times from a state i can be simply defined 
as the probability to stay for n — 1 steps in state i times the probability to 
exit at step n, i.e.: 

P l {n) = {l-W ll ){W ll ) n - 1 . (33) 
that can be rewritten as: 



Pi(n) 



exp(— na) , with a = | In Wu\ . (34) 



In a similar way we can define the probability fijin) of the first arrival 
from state i to state j at step n. This is the probability to arrive to state 
j starting from i in n step i.e. (W n )ij minus the probability of first arrival 
at step n — k times the probability of return in k steps, i.e. {W k )jj with 
k — 1, • • ■ , n — 1: 



n-1 



/ y (n) = (W^- - E - k)(W%. (35) 



fe=i 



For each state of a Markov process a correlation function can be defined 
for the variable Xi{ n ) which is equal to 1 if at time n the i state is visited 
and to zero otherwise (see eqs. |23| - p4D . The normalized correlation function 



= < xM) > - < xM >\ (36) 

is strictly related to the diagonal element (W n )u of the transition matrix to 
the n-th power and to the stationary frequency Pi. Notice that 

< X *(0)>= P i7 < Xl (0) 2 >= P. (37) 

Furthermore, being Pj the probability that the initial state at n = be i 
and (W n )a the probability to be in i again after n iterations, one has: 

<Xi{0)Xi(n) >= P t -(W n ) u (38) 

and therefore: 

Q(n) = ((W n ) tl -P)/(l-P)) (39) 
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TABLES 



Table 1 

Transition matrix elements 





Case A 


Case B 


Case C 




.66 


.74 


.58 


w 12 


.17 


.13 


.21 


w 13 


.17 


.13 


.21 


w 21 


.12 


.09 


.14 


W22 


.88 


.91 


.86 


w 23 


.00 


.00 


.00 


w 31 


.12 


.09 


.14 


w 32 


.00 


.00 


.00 


w 33 


.88 


.91 


.86 



Table 2 

Visit probabilities 



Pi 


Case A 


Case B 


Case C 


Pi 


.26 


.25 


.25 


P2 


.37 


.375 


.375 


P3 


.37 


.375 


.375 



Case A - Deterministic chaotic model defined by eq. (^) related to the 
stream function (0) with parameters L = 7.5, B = 1.2, c = 0.12, uj = 
0.4, e = 0.3. 

Case B - Turbulent diffusion model defined by eqs. ( pT| - [2"8| ) with param- 
eters o = .05 t = T/4. 

Case C - Model with chaotic advection plus turbulent diffusion with the 
same parameters of case A and B. 

The statistics have been computed over 2 10 6 periods. 
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FIGURE CAPTIONS 



FIGURE 1: Snapshot of the velocity field derived from the stream 
function eq.(|) with L = 7.5, B = 1.2, c = 0.12. 

FIGURE 2: Stream lines of the time-dependent stream function eq.(|2]), 
with B given by eq.(|), B = 1.2, to = 0.4 and e = 0.3 (T = 2it/lu), at 
three different times: (a) t = 0. (T/2), (b) t = T/4 and (c) t = 3T/4. 

FIGURE 3: Critical values of the periodic perturbation amplitude for 
the overlap of the resonances, € c /Bq vs uj/ujq, for the stream function 
d) with L = 7.5, B = 1.2, c = 0.12 and cu = .25, which is the 
typcal frequency for the rotation of a tracer on the boundary of the 
recirculation gyres. The critical values have been estimated following 
a cloud of 100 particles initially located between the states 1 and 2 up 
to 500 periods. 

FIGURE 4: Spreading of a cloud of 5000 tracers at different times for 
the deterministic model (see Fig. 2): (a) t — all tracers are inside a 
very small square, (b) t = T ,(c) t = 5T, (d) t = 10 T, (e) t = 20 T 
and (e) t = 100 T. 

FIGURE 5: Probability distribution of the first exit times from states 
1 (diamonds), 2 (squares) and 3 (crosses) for the deterministic model 
(see Fig. 2). The straight lines are the markovian predictions given by 
eq.(^) with Wa of table 1 (case A). The time unit is the period T of 
the perturbation. The statistics is computed over 2 10 6 periods. 

FIGURE 6: Correlation functions of the states 1 (diamonds) and 2 
(crosses) compared with the markovian predictions (continuous lines) 
eq. (p9|) for the deterministic model (see Fig. 2). 

FIGURE 7: Probability distribution of the meridonal mixing times 
(MMT) compared with the markovian predictions (continuous lines) 
eq. (|3~3l) for the deterministic model (see Fig. 2). 

FIGURE 8: Comparison beetwen the symbolic sequence of the states 
as function of time (for t = T, 2T, • • •) obtained by the integration of 
the deterministic model equations (1,2) (dotted lines) and the symbolic 
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sequence generated from the Markov chain (solid lines) defined by the 
transition matrix computed as described in eq.(|i~8D and reported in 
table 1 (case A) . Here represents the state 1 i.e. the jet while 1, — 1 
the recirculation gyres i.e. the states 2,3. 



FIGURE 9: Block entropies h n vs n flT5| ) for the deterministic model 
(see Fig. 2), computed from a sequence of 10 6 symbols. 

FIGURE 10: The same as figure 4 for the stochastic model given by 
eqs. (|27| , |28|) , with parameters a = 0.05 and r = T/4. The time unit T 
is set equal to the period of the deterministic perturbation (see Fig. 2). 

FIGURE 11: Probability distribution of the first exit times from states 
1,2 and 3 for the stochastic model (see Fig. 10). The straight lines are 
the markovian predictions given by eq.fl3"3]) with Wu of table 1 (case 
B). 

FIGURE 12: Correlation functions compared with the markovian pre- 
dictions (continuous lines) eq. (|39|) for the stochastic model (see Fig. 
10). 

FIGURE 13: Probability distribution of the MMT compared with the 
markovian predictions eq. (|3~5D for the stochastic model (see Fig. 10). 

FIGURE 14: Block entropies h n vs n for the stochastic model, com- 
puted from a sequence of 10 6 symbols. 

FIGURE 15: Probability distributions of the first exit times from states 
1,2 and 3 in the model with chaotic advection combined with turbulent 
diffusion (sect. 4.3) with parameters: B = 1.2, u = 0.4, e = 0.3 
and a = 0.05, r = T/4 where T = 2-k/uj. The straight lines are the 
markovian predictions given by eq. ([33]) with Wu of table 1 (case C). 



FIGURE 16: Probability distributions of the MMT compared with 



the markovian predictions eq.(35) for the model and the parameters of 
figure 15. 

FIGURE 17: Correlation functions compared with the markovian pre- 
dictions eq. (|39|) for the model and the parameters of figure 15. 
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.919E+00 

> 

MAXIMUM VECTOR 




CONTOUR FROM -.90000 TO .90000 CONTOUR INTERVAL OF .10000 PT(3,3)= .60931 




CONTOUR FROM -.80000 TO .80000 CONTOUR INTERVAL OF .10000 PT(3.3)= .60922 




CONTOUR FROM -.80000 TO .80000 CONTOUR INTERVAL OF .10000 PT(3.3)= .60906 
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